Order parameters in the Verwey phase transition 
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The Verwey phase transition in magnetite is analyzed on the basis of the Landau theory. The free 
energy functional is expanded in a series of components belonging to the primary and secondary 
order parameters. A low-temperature phase with the monoclinic P2/c symmetry is a result of 
condensation of two order parameters X3 and A5. The temperature dependence of the shear elastic 
constant C44 is derived and the mechanism of its softening is discussed. 
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I. INTRODUCTION 

The mechanism of the phase transition in magnetite (Fe304) at Ty = 122 K, discovered by Verwey has 
remained a big puzzle in the condensed matter physics for almost 70 years. Developments in experimental and 
theoretical methods during last years enabled to reveal subtle changes in the crystal and electronic structure below 
Ty @- S] • A simple charge ordering picture in which metal-insulator transition is induced by electrostatic interactions 
was replaced by a highly complex scenario in which lattice, charge, spin and orbital degrees of freedom are involved. 
Recent theoretical studies revealed the important role of the local electron interactions and orbital correlations in the 
t2 g states on iron ions Q. 

The electronic interactions are complemented by the lattice deformation, which breaks the cubic symmetry and 
induces a low-temperature (LT) monoclinic phase driven by the electron-phonon interactions Q ■ In the previous work 
[El IB] j w e have obtained the phonon spectrum of magnetite (Fig. 1) using the ab initio computational technique 0]. 
We have identified two primary order parameters (OPs) at kjc = ^(0, 0, 1) and kA = ^-(0, 0, \) with the X 3 and A 5 
symmetry, respectively, which both play important role in the VT: (i) the A5 mode is responsible for the doubling 
of the unit cell along the c direction in the monoclinic phase, while (ii) the A3 phonon induces the metal-insulator 
transition by its coupling to the electronic states near the Fermi energy [5[. Due to the electron-phonon interaction 
the above OPs are combinations of the electron (charge-orbital) and lattice components. This explains why the 
phonon soft mode has not been observed. Instead, low-energy critical fluctuations of OPs were found by the diffuse 
neutron scattering . The condensation of the OPs below Ty explains the crystal symmetry change as well as the 
charge-orbital ordering. 

The group theory predicts also secondary OPs, which do not effect the symmetry below Ty but modify the properties 
of magnetite close to a transition point. At the T point, the T 2g mode can be classified as the secondary OP, and its 
coupling to the shear strain explains the softening of the C44 elastic constant [101. The lowest T 2g optic mode, 
marked in Fig. 1, could contribute quantitatively to the free energy, but it does not play any significant role for the 
VT. 

In this work, we introduce and analyze the Landau free energy for the VT, and discuss a solution corresponding to 
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FIG. 1: The phonon dispersion curves of Fe304 (solid lines) compared with the experimental points taken from Ref. Q. The 
order parameters X3, A5, and Ti g are marked by circles. 
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the LT monoclinic phase. We derive also the temperature dependence of C44. 



II. FREE ENERGY 



The Landau free energy can be expanded into a series of the components of the OPs. The invariant terms describing 
couplings between the OPs were derived using the group theory methods [T2I ]. The only nonzero components of the 
primary OPs X 3 and A 5 are denoted by g and q, respectively. We include also the secondary OP with the T 2g 
symmetry (77) and shear-strain (e). The free energy can be written in the form [|| 

T -r- , "1 2 , Pi 4 , 71 6 , "2 2 , #2 4 , *1 2 2 , a 3 2 , a i 2 , S 2 2 S 3 2 

^ = J r o+y9 + —g +-q 9 + ~2 Q ~ q ~2 9 Q Y V Y € ' + ~2 VQ Y £g +<54??e ' W 
were Tq is the part of the potential, which does not change through the transition. We assume that j3\ > 0, 02 > and 
7i > to ensure the stability of the potential at high temperatures. For the second-order terms we assume standard 
temperature behavior ai = a,i(T — T C i) near the critical temperature T C i for i — 1, 2, 3, which would correspond to a 
continuous phase transition. The coefficient 0J4 is the shear elastic constant at high temperatures (C44). The coupling 
between the primary OPs is biquadratic, between the secondary and primary OPs has the linear-quadratic form, and 
the coupling between the components of the secondary OP is of the bilinear type. Taking first derivatives of T over 
all OPs we get 



g(a x + (3 ig 2 + 7l .g 4 + S iq 2 + 5 2 r) + S 3 e) = 0, (2) 
q(a 2 + foq 2 + 6 l9 2 ) = 0, (3) 



dg 
dT_ 
dq 

d S 

— = a 3 r/ + 6 t e + y g 2 = 0, (4) 

— = a 4 e + S^ + j-g 2 = Q. (5) 

The solution g = q = ri = e = corresponds to the high-temperature cubic symmetry (Fd3m). From Eq. (3) we 
obtain the dependence between g and q 

2 S ig 2 + a 2 . , 

1 = 3 > (6) 

which has three possible solutions: (i) g = and q 2 = — |g if a 2 < (Pbcm), (ii) q = and g 2 = — if a 2 > 
and Si > or a 2 < and Si > (Pmna), (Hi) g ^ and 5 / (P2/c). In the brackets we put the space group 
symbols, which characterize the low-symmetry phases. The solution (Hi) which corresponds to the experimentally 
observed LT monoclinic phase requires simultaneous condensation of both primary OPs. The necessary condition for 
this is a negative value of S±. Indeed, it has been established by the ab initio studies that the total energy is lowered 
when the crystal is distorted by both X 3 and A5 modes |f>j. For Si < 0, Eq. ([6]) has a non-zero solution provided that 
\S1\g 2 > a 2 . It implies that for a 2 > (T > T c2 ), the phase transition occurs when the OP g exceeds a critical value 
t^t, so it has a discontinuous (first-order) character. 

From Eqs. (4) and (5) we get 

S 3 S 4 - S 2 a 4 2 \ 2 <*2&i - S 3 a 3 2 2 

2a 3 a<i — 2a\ 2a 3 a^ — 2<5| 

which shows that t] ^ and e ^ only if <? 7^ 0. Eliminating g, rj and e using Eqs. (6) and (7), the potential T can 
be written as a function of g 

T T> 1 a 2 , P 4 , 71 6 /o\ 

•^ = •^0 + ^.9 +^9 , (8) 

where the renormalized coefficients are 

H=Fo~ J"> (9) 
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a = ai — , (10) 
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where 



The zero-order and second-order terms depend on the parameters belonging to the primary OPs. The secondary OPs 
modify only the forth-order term. In this notation, the solution of Eqs. (2)-(5), which minimizes the potential T 
reads 

27 Pi 

To study the softening of C44, we have eliminated g, q, and rj using Eqs. (2), (3), and (4), and expressed the free 
energy as a function of e only. In these calculations we have omitted the sixth-order term, which usually has a small 
contribution near the transition point. The elastic constant C44 is obtained using the standard definition 

d 2 T S 2 5 2 

C44(T)^ — = C 4 °4-^-^-, (13) 

5 = 5,-^, a > 3 = a 3 -^-=a 3 (T-T> 3 ), = ft - (14) 

2pi Wi lh 

with T' c3 = T c3 + S 2 /2(3[. The second and third term in Eq. (fT3| are negative at high temperatures, so both 
contribute to the softening of C44. It means that all couplings included in Eq. |T|) are involved in this behavior. 
The main temperature dep endence is caused by the second term, but also the last term in Eq. Q13[) may depend on 
temperature. In Ref. 10], the softening of C44 was explained taking into account only the bilinear coupling (64). 
Indeed, if we assume that 62 = 5 3 = 0, Eq. (fl~3|) reduces to that found in Ref. C44 = C44 — S 2 /a 3 . Omitting the 
last term, Eq. (|13|) can be written in the form 

C44 = C*4 4 - - T ° , (15) 

where To = T' ci + S 2 / 6440,3. Since, Eq. (fTS)) has the same form as that one discussed in Ref. [l(|, the fitting procedure 
will give the following values of parameters: To = 66 K and T' c3 = 56 K. Note that the meaning of both temperatures 
T and T' c3 is here different than in Ref. [l(| as they include all interactions. The softening is not complete (only 
about 10 %) since the first-order phase transition and the change of structure occurs at much higher temperature 
than Tq. 



III. DISCUSSION 



The present work shows that the VT can be analyzed using the Landau theory of phase transitions. It provides a 
basis to study the interplay between the OPs and the mechanism of phase transition with two or more OPs. Moreover, 
majority of experimental facts can be understood in one coherent picture. We emphasize that the coupling between 
the modes with the X3 and A5 symmetry plays the crucial role in the VT in magnetite and explains both the 
occurrence of the LT monoclinic phase and the metal-insulator transition. This results in rather complex structure 
of the charge and orbital ordering, with two different charge modulations characterized by kx and kA wave vectors 
found in diffraction studies 2]. In this work, we have discussed the conditions, which have to be fulfilled for both 
OPs to develop simultaneously. It depends primarily on the coupling coefficient Si, which should be negative below 
Ty and large enough to stabilize both OPs. 

Our present analysis extends also the previous studies of the temperature dependence of C44. So far, the softening 
of C44 was studied in the models restricted to the T point [l(| EH • Since the phase transition is driven mainly by the 
OPs at k / 0, the OPs at the zone center develop solely due to their coupling to other modes. Therefore, the bilinear 
coupling at the T point is rather a side effect of the phase transition, not the main origin of it. The present study 
shows that also a direct coupling between the X3 primary OP and the shear strain (S3) influences the C44 elastic 
constant, providing an alternative mechanism of its softening. 
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